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Time-dependent density functionals in principle depend 
on the initial state of a system, but this is ignored in func- 
tional approximations presently in use. For one electron, it 
is shown that there is no initial-state dependence: for any 
density, only one initial state produces a well-behaved poten- 
tial. For two non-interacting electrons with the same spin in 
one-dimension, an initial potential that makes an alternative 
initial wavefunction evolve with the same density and current 
as a ground state is calculated. This potential is well-behaved, 
and can be made arbitrarily different from the original poten- 
tial. 



I. INTRODUCTION AND CONCLUSIONS 



Ground-state density functional theory |l|,|2j has had 
an enormous impact on solid-state physics since its in- 
vention, and on quantum chemistry in recent years ||. 
Time-dependent density functional theory (TDDFT) al- 
lows the external potential acting on the electrons to be 
time-dependent, and so opens the door to a wealth of 
interesting and important phenomena that are not eas- 
ily accessible, if at all, within the static theory. Impor- 
tant examples include atomic and molecular collisions 
jjj, atoms and molecules in intense laser fields ||, elec- 
tronic transition energies and oscillator strengths HQ, 
frequency-dependent polarizabilities and hyperpolariz- 
abilities, etc. ||, and there has been an explosion of time- 
dependent Kohn-Sham calculations in all these fields. 
In almost all these calculations, the ubiquitous adia- 
batic local density approximation (ALDA) ||[l0| is used 
to approximate the unknown time-dependent exchange- 



correlation potential, i.e., v 



ALDAr 



i](rt) 



,,unif 



(n(rt)), 



where i>™ (n) is the ground-state exchange-correlation 
potential of a uniform electron gas of density n. While 
this seems adequate for many purposes [0, little is 
known about its accuracy under the myriad of circum- 
stances in which it has been applied. 

Runge and Gross |L2| formally established time- 
dependent density functional theory (TDDFT), show- 
ing that for a given initial state, the evolving den- 
sity uniquely identifies the (time-dependent) potential. 
This established the correspondence of a unique non- 
interacting system to each interacting system and so a 
set of one-particle Kohn-Sham equations, much like in 
the static theory. This one-to-one mapping between den- 
sities and potentials is the time-dependent analog of the 
Hohenberg-Kohn theorem, but with a major difference: 



in the time-dependent case, the mapping is unique only 
for a specified initial state. The functionals in TDDFT 
depend not only on the time-dependent density but also 
on the initial state. This dependence is largely unex- 
plored and indeed often neglected, for example in the 
ALDA for the exchange-correlation potential mentioned 
above. 

What do we mean by an initial-state dependence? In 
the ground-state theory, there is a simple one-to-one re- 
lation between ground-state densities and Kohn-Sham 
potentials i? s , assuming they exist. For example, for 
one electron in one dimension, we can easily invert the 
Schrodinger equation, to yield 



v s (x) 



d 2 y / n(x)/dx 2 
2^/nJx) 



(1) 



where n(x) is the ground-state density. We use atomic 
units throughout (Ti = m — e 2 = 1). For N elec- 
trons in three dimensions, one can easily imagine contin- 
uously altering v s (r), the Kohn-Sham potential, solving 
the Schrodinger equation, finding the orbitals and cal- 
culating their density, until the correct v s (r) is found 
to reproduce the desired density. By the Hohcnbcrg- 
Kohn theorem, this potential is unique, and several clever 
schemes for implementing this idea appear in the litera- 
ture |l3[|l5j -|2l[] . This procedure could in principle be im- 
plemented for interacting electrons, if a sufficiently versa- 
tile and accurate interacting Schrodinger equation-solver 
were available. 

Now consider the one-dimensional one-electron density 



n(x) = 2x 2 exp(— a: 2 )/i/7r, 



(2) 



(actually the density of the first excited state of a har- 
monic oscillator) . If we consider this as a ground-state 
density, we are in for an unpleasant surprise. Feeding 
it into Eq. ([!]), we find that the potential which gener- 
ates this density is parabolic almost everywhere (a; 2 /2), 
but has a nasty unphysical spike at x = 0, of the form 
(5(x)/|x|. We usually exclude such potentials from con- 
sideration p2[ , and regard this density as not being v- 
representable. 

But now imagine that density as being the density of a 
first excited state. In that case, the relation between den- 
sity and potential is different, because the orbital changes 
sign at the node. The mapping becomes 
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v s (x) = 



(sgn(x - x )^/n(x) \ /da 



2 (sgn(a; - x Q ) v / n(x) 



(3) 



where sgn(a;) = 1 for x > and -1 for x < 0, and n(xo) = 
0. If we use this mapping, we find a perfectly smooth 
parabolic well (x 2 /2). This is a simple example of how 
the mapping between densities and potentials depends on 
the initial state. 

More generally, for any given time-dependent density 
n(rt), we ask how the potential v(rt) whose wavefunc- 
tion yields that density depends on the choice of initial 
wavefunction ^q, i- e -j m general v^o, n] (rt). Our aim in 
this paper is to explicitly calculate two different poten- 
tials giving rise to the same time-dependent density by 
having two different initial states. Note that even finding 
such a case is non-trivial. The choice of wavefunctions 
is greatly restricted by the time-dependent density. As 
van Leeuwen has pointed out (23), the continuity equa- 
tion h = — V ■ j implies that only wavefunctions that 
have the correct initial current are candidates for gener- 
ating a given time-dependent density. Van Leeuwen also 
shows how to explicitly construct the potential generat- 
ing a given density from an allowed initial wavefunction 
using equations of motion. 

Why is this important? The exchange-correlation po- 
tential, v KC (rt), of TDDFT is the difference between a 
Kohn-Sham potential and the sum of the external and 
Hartree potentials. Since both the interacting and non- 
interacting mappings can depend on the choice of initial 
state, this potential is a functional of both initial states 
and the density, i.e., v xc [n, ^q, $o]( r O- But in common 
practice, only the dependence on the density is approxi- 
mated. We show below that this misses significant depen- 
dences on the initial state, which can in turn be related 
to memory effects, i.e., dependences on the density at 
prior times. 

In the special case of one electron, we prove in sec- 
tion II A that only one initial state has a physically well- 
behaved potential. Any attempt to find another initial 
state which evolves in a different potential with the same 
evolving density, results in a "pathological" potential. 
The potential either has the strong features at nodes 
mentioned above, or rapidly plunges to minus infinity 
at large distances where the density decays. (How such 
a potential can support a localized density is discussed 
.) Such non-physical states and potentials 



in section III 



are excluded from consideration (as indeed they are in 
the Runge-Gross theorem) . Thus there is no initial-state 
dependence for one electron. 

We might then reasonably ask, can we ever find a well- 
behaved potential for more than one allowed initial wave- 
function? The answer is yes, which we demonstrate with 
a specific example. Consider two non-interacting elec- 
trons of the same spin in a harmonic well. In the ground 
state of this two electron system, the first electron occu- 




FIG. 1. The top left-hand plot shows the ground-state or- 
bitals <f>i (solid) and 02 (dashed) and their density n (thick 
solid line) for the harmonic potential in the lower left-hand 
plot. The top right-hand plot contains the real and imagi- 
nary parts of alternative orbitals 0i (solid) and 02 (dashed), 
and their density, n(thick solid), while below is the unique 
initial potential v (solid line) that keeps the density constant. 



pies the oscillator ground state, and the second occupies 
the first excited state, as shown in Fig. [|. If we keep 
the potential constant, the density will not change. By 
multiplying each orbital by a spatially varying phase, and 
choosing these phases to make the current vanish, we find 
an allowed alternative initial state (see section |n] for de- 
tails). Van Leeuwen's prescription then yields the unique 
potential which makes this wavefunction evolve with the 
same density. The difference is perfectly well-behaved, 
and can be made arbitrarily large by adjusting a con- 
stant in the phases of the alternative orbitals. To our 
knowledge, this is the first explicit construction of two 
different potentials that yield the same time-dependent 
density. Other examples are given in section [II. 

Now imagine the density of Fig. | were the ground- 
state density of some interacting two-electron system, in 
some external potential v cx t(x)- Then both potentials 
shown in lower panels of Fig. [I] are possible Kohn-Sham 
potentials, v s (x), for this system. Since the Hartree po- 
tential is uniquely determined by the density, we have 
two very different exchange-correlation potentials, differ- 
ing by the amount shown. In fact, different choices of 
initial wavefunction allow us to make the two dips arbi- 
trarily deep or small. Any purely density functional ap- 
proximation misses this effect entirely, and will produce 
the same exchange-correlation potential for all cases. 

So, even in the simplest case of non-degenerate inter- 
acting and Kohn-Sham ground states, one can choose 
an alternative Kohn-Sham initial state, whose potential 
will look very different from that which evolves from the 
initial ground state. In practice, the majority of applica- 
tions of TDDFT presently involve response properties of 
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the ground-state of a system, and one naturally choses 
to start the Kohn-Sham system in its ground state. This 
choice is also dictated by the common use of adiabatic ap- 
proximations for exchange-correlation potentials, which 
are approximate ground-state potentials evaluated on the 
instantaneous density. Such models will clearly be inac- 
curate even at t = if we start our Kohn-Sham calcula- 
tion in any state other than its ground state. 

The initial-state dependence of functionals is deeply 
connected to the issue of memory effects which are ig- 
nored in most TDDFT functional approximations used 
today. Yet these can often play a large role in exchange- 
correlation energies in fully time-dependent (i.e. non- 
perturbative) calculations [53 as well as giving rise to 
frequency-dependence of the exchange-correlation kernel 
in linear response theory (f xc (ui)) |2j|. Functionals in 
general depend not only on the density at the present 
time, but also on its history. They may have a very non- 
local (in time) dependence on the density. But still more 
about the past is required: the functional is also haunted 
by the initial wavefunction. The initial state dependence 
is inextricably linked to the history of the density and 
in fact can be absorbed into density-dependence along a 
pseudo pre- history Q . The results of this current pa- 
per shed some light on the importance of memory effects 
arising from the initial wavefunction. 

To summarize, we have shown that there is no initial- 
state dependence for one electron, and that there can 
be arbitrarily large initial-state dependence for two elec- 
trons. 



II. THEORY 

Consider a density n(rt) evolving in time under a 
(time-dependent) potential v(rt). Can we obtain the 
same evolving density n(rt) by propagating some differ- 
ent initial state in a different potential v(rt)7 This was 
answered in the affirmative in Ref. |23| under the condi- 
tion that the two initial states have the same initial den- 
sity and initial first time-derivative of the density. Here 
we shall show that additional restrictions are required on 
the initial state for this statement to hold. In the one- 
electron case, the additional restrictions are so strong 
that there is no other initial state that evolves with the 
same density as another does in a different potential. 



A. One electron 

Any two one-electron wavefunctions <p(rt) and <j>(rt) 
with the same density n(rt) = \<fi(rt)\ 2 = \4>(rt)\ 2 , are 
related by a spatial and time dependent phase factor: 



where the phase a(rt) is real. The evolution of 
each wavefunction is determined by the time-dependent 
Schrodinger equation with its potential (dot implies a 
time derivative): 

(-V 2 /2 + v(rt)) (f>(rt) = i<p(rt). (5) 

Both will satisfy the continuity equation: 

n(rt) = -V-j(rt) (6) 

where the current density of a wavefunction <fi is 

j(rt) = i fa(rt)V0*(rt) - 0*(rt)V#rt)) /2. (7) 

Substituting 4>(rt) from Eq. (|) into Eq. we obtain 

Aj(rf) = n(rt)Va(vt). (8) 

(We use the notation Aa to denote a — a.) Because the 
densities are the same for all times, An(rt) = 0, so by 
Eq. © 



V • [n(rt)Va(rt)] = 



(9) 



Integrating Eq. (^J) with a(rt) and performing the inte- 
gral by parts, we find 



J d 3 r n{rt)\\7a(rt)\ 2 = 0. 



(10) 



4>(rt) = 4>(rt) exp (ia(rt)) 



(4) 



We have taken the surface term J d 2 S ■ (cmVa) evalu- 
ated on a closed surface at infinity, to be zero: this arises 
from the physical requirement that at infinity, where the 
electron density decays, any physical potential remains 
finite. (In fact this condition is required in the proof 
of the Runge-Gross theorem |12|). If the surface term 
did not vanish, then aVa must grow at least as fast as 
(r 2 n(r)) -1 as r approaches infinity. This would lead to a 
potential that slides down to — oo which can be seen by 
inversion of the time-dependent Schrodinger equation in 
the limit of large distances. The state would oscillate in- 
finitely wildly at large distances in the tails of the density 
but the decay of the density is not enough to compen- 
sate for the energy that the wild oscillations impart: this 
state would have infinite kinetic energy, momentum and 
poten tial energy. (We shall see this explicitly in section 
[II A ). So for physical situations, the surface term van- 
ishes. 

Because the integrand above cannot be negative, yet 
it integrates to zero, the integrand itself must be zero ev- 
erywhere. Thus Va(rt) = everywhere except perhaps 
at nodes of the wavefunction where n(roi) = 0. In fact, 
even at the nodes, Va(roi) = to avoid highly singular- 
potentials: if Vq was finite at the nodes and zero every- 
where else, then as a distribution, it is equivalent to being 
zero everywhere, for example its integral a is constant. 
There remains the possibility that Vq is a sum of delta 
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functions centered at the nodes, however this leads to po- 
tentials which are highly singular at the nodes, as in the 
introduction. Such unphysical potentials are excluded 
from consideration, so that Va(rt) = 0, or a(rt) = c(t). 
The wavefunctions <p(rt) and <p(rt) can therefore differ 
only by an irrelevant time-dependent phase. In particu- 
lar, this means that only one initial state and one poten- 
tial can give rise to a particular density, i.e, the evolving 
density is enough to completely determine the potential 
and initial state. 

The one-electron case is a simple counter-example to 
the conclusions in Ref. which rely on the existence 
of a solution to 



V • [nVAu] = r](rt), 



(11) 



where rj(rt) involves expectation values of derivatives of 
the momentum-stress tensor and derivatives of the in- 



teraction (see section [IB ). This is to be solved for the 
potential subject to the requirements that the two initial 
states have the same n(rQ) and n(r0), and that VAv — > 
as r — > oo. The two initial wavefunctions in the one- 
electron case have the same initial n(rO) (Eq. [|) and 
n(rO) (Eq. ^) but no two physical potentials exist under 
which they would evolve with the same density, because 
there is no solution to Eq. ([ll]) subject to the boundary 
condition that VAu = 0. (For an explicit demonstra- 
tion of this in one-dimension, see section III A.) We shall 



come back to this point at the end of section 



II B 



B. The many-electron case 

In this section we follow van Leeuwen's prescription to 
find the potential needed to make a given initial state 
evolve with the same density as that of another. How- 
ever we simplify the equations there somewhat to make 
the search for the solution of the potential easier. Given 
an initial state \&o which evolves with density n(rt) in a 
potential v(rt), we solve for the potential v(rt) in which 
a state *f> evolves with the same density n(rt). If we re- 
quire \1/ to have the same initial density and initial first 
time-derivative of the density, then a solution for v may 
be obtained from equating the equations for the second 
derivatives of the density for each wavefunction, subject 
to an appropriate boundary condition like Av — > at 
large distances. We are not guaranteed that such a so- 
lution exists: the wavefunction must have the additional 
restriction that the initial potential computed in this way 
is bounded at infinity. 

The equation of motion for h yields (Eq. (15) of Ref. 

PI): 



ii{rt) = V • [n(rt)Vv(rt) + t(ri) + f ee (ri)] (12) 



where 



t(rf) - (V - V)(V 2 - V 2 ) Pl (r'rt)\ r , =r /4 (13) 

and 

fee(rt) - J d 3 r'P(v'rt)Vv cc (\r' - r|)/2 (14) 

where pi(r'rt) is the (off-diagonal) one-electron reduced 
density matrix, P(r'ri) is the pair density (diagonal two- 
electron reduced density matrix) and v oc (u) is the two- 
particle interaction, e.g. 1/u. Here and following, V and 
V indicate the partial gradient operators with respect to 
r and r' respectively. In Eq. ([l3]) and similar following 
equations, r' is set equal to r after the derivatives are 
taken. 

The idea J2i| is to subtract Eq. ( |l2| ) for wavefunction 
^ from that for wavefunction "J , and require that n is the 
same for each. First, we simplify the kinetic- type term, t. 
Differentiating the continuity equation (Eq. (^)) implies 



(V + V')(V 2 - V' 2 )A Pl (r'rt)| r ^ r = 



(15) 



This equality enables us to incorporate the satisfaction of 
the equation of continuity in Eq. ([l2]) (when we subtract 
the equation for 'J from that for and it also simplifies 
the kinetic-type term: 



At(rt) = -V (V 2 - V' 2 ) Ap 1 (r'rt)| P /= r /2. 



(16) 



(We note that although this is no longer explicitly real, 
it is in fact real for states with the same density and 
first time-derivative.) So our simplified equation to solve 
becomes 



V • [nVAu + At + Af cc ] = 



(17) 



where At is given by Eq. (|l6|) and Af ee is given by Eq. 
fll4|), applied to the pair density difference. 

To calculate the derivatives in the kinetic-type term, 
we define 



7(rr'i) = A log p\ (r'rt) 



(18) 



Note that 7 vanishes at r = r', and since pi(r'ri) = 
p*(rr't), 7(r'rf) = 7*(rr'i). These relations also imply 
that V m 7(r'rt)| r , =r = V m 7*(rr'i)| r , =r . Writing 



7(rr'i) = /?(rr'i) + ia(rr'i) 



(19) 



where a and /3 are real functions, we also find 
V/3(r'rt)| r / =r = 0, since V r [/3(r'ri)| r > =r ] = 0. Also, 
V • V'a(r'rt)| r / =r = 0, which follows from the antisym- 
metry of a. The generalization of Eq. (0) is 



Aj(rt) = rc(rt)Va(r'rt)| r / =r 



(20) 



Continuity (Eq.(ra)) then gives us a condition on the near- 
diagonal elements of a: 



V ■ [n(ri)Va(r'ri)| P , =r ] = 



(21) 
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Using all these results in Eq. (|l7|), we find 



Av = ia"/2 + ia'<f>'/(f> - a' /2 - a = 0. (23) 



VAv(rt) = (nVa - V 2 aV3pi + V x B) /n 
+ V3? (V 7 • Vpi - V'7 • V'pi) In 

+ \™ ((V7) 2 - (V'7) 2 ) + \v (V 2 - V' 2 ) (3 
~ 2nJ^T) J d 3 r'AP(rr't)Vv cc (W ~ r|) (22) 

where in the first three lines we have omitted the argu- 
ments and it is understood that r' and r are set equal 
after all the derivatives are taken. B(rt) is an unde- 
termined vector whose role, together with an additional 
constant C(t), is to ensure satisfaction of a boundary 
condition on the potential. 

Now the prescription is to pick an initial state which 
has the same initial density and initial first time- 
derivative of the density as the state "J; that is, require 
7(rr0) = and Eq. taken at t = 0. Then one can 
evaluate Eq. ( |22] ) at t = and so find v(r0). The pro- 
cedure for t > is described in detail in Ref. ^3|. In 
order for this procedure to yield a well-behaved physical 
potential, one needs to first check that the initial poten- 
tial is not divergent at infinity. (Equivalently, we may 
require that the elements of the momentum-stress ten- 
sor appearing in Eq. ([ll]) do not diverge at infinity). 
This gives an additional restriction on the initial state. 
In the one-electron case, this restriction rules out any 
other candidate for an initial wavefunction which evolves 
with the same density as another wavefunction does in 
another potential: there is no way to pick B(r0) or the 
constant C(0) to satisfy any physical boundary condi- 
tion discussed above. In the many-electron case, our ad- 
ditional condition restricts the allowable wavefunctions, 
but does not render the question of initial-state depen- 
dence moot as in the one-electron case. 



III. EXAMPLES 

A. One electron in one dimension 

By studying the time-dependent Schrodinger equation 
for one electron in one dimension, it is simple to find ex- 
plicitly the potential v(xt) which cajoles <f>(xt) into evolv- 
ing with the same density as 4>(xt) which evolves in a 
different potential v(xt). Consistent with the conclusions 
above, this potential diverges to —00 at large x, which is 
unphysical. The initial state is pathological in the sense 
that its expectation value of momentum, kinetic energy 
and potential energy all diverge. A phase-space picture 
helps us to see how such a potential can hold a localized 
density. 

Inserting Eq. (0) into the time-dependent Schrodinger 
Eq. (0) and calculating the derivatives, we obtain 



where primes denote spatial derivatives. We now write 
in terms of an amplitude and phase 



(xt) = \f n(xt) exp(i9(xt)). 



(24) 



Substituting into Eq. (|23| ) and setting the real and imag- 
inary terms separately to zero, yields 

Av = -a-a'9' - a' 2 /2; a" + a'n'/n = 0. (25) 



We find for a': 



c{t)/n{xt). 



(26) 



We observe that this is also obtained when Eq. (g) (which 
arose from setting the time-derivatives of the densities 
to be equal) is considered in one-dimension. Integrating 
once more gives 



i(xt) = c(t) 



dx' 



n(x't) 



Plugging this solution into Eq. 



d(t) 



gives 



(27) 



Av(xt) 



c{t)0'{xt) 1 c 2 {t) 



-m 



n{xt) 2n 2 (xt) 

dx' , , f x n(x't) , , ;. . . . 



n(x't) 



n 2 (x't) 



We see immediately the divergence of this potential 
(for non-zero e(t)) where the density decays at large 
x. This demonstrates by explicit solution of the time- 
dependent Schrodinger equation that the only potentials 
in which a density can be made to evolve as the density 
in some other potential a re u nphysical, consistent with 
the conclusions of section [I A . 

The state <fi oscillates more and more wildly as x gets 
larger in the tails of the wavefunction. Although the de- 
cay of the density at large distances unweights the rapidly 
oscillating phase, it is not enough to cancel the infinite 
energy that the wild oscillations contribute. Calculating 
the expectation value of momentum or kinetic energy in 
the state Eq. (||) with a given by Eq. (p^j), for a typical 
density and state <f> (e.g., one which decays exponentially 
at large x), we find they blow up. 

It may be at first glance striking that a potential which 
plunges to minus infinity at large distances can hold a 
wavefunction which is localized in a finite region in space. 
Consider the special case in that 4>(xi) is an eigenstate 
of a time-independent potential v(x). Let us also choose 
c(t) = c to be time- independent, so that v(x) is also 
time-independent and <p(xt) is an eigenstate of it. Let 
the density n(x) = \<fi(xt)\ 2 = \4>(xt)\ 2 be localized at the 
origin. For example, v might be a potential well with 
flat asymptotes. Then we have the interesting situation 
where the eigenstate 4>(x) is localized at the origin of its 
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FIG. 2. The lower figure shows classical phase-space con- 
tours for the sech 2 well. The top figure shows the potential 
(dashed line), the wavefunction (solid line) corresponding to 
the heavily drawn contour in the phase-space below and its 
density (thick solid line). 

potential v which plummets to — oo at large x. In Figs. 
H and H we have plotted the potentials, the densities and 
the (real part of the) wavefunctions for the two cases; 
notice the steep cliffs of v and the rapid oscillations of </> 
as x gets large as we predicted. We chose the potential 
v = — sech 2 (x), and the state <f> as its ground-state. In 
the lower half of each figure are the classical phase space 
pictures (the classical energy contours) for the two poten- 
tials; to a good approximation the quantum eigenstates 
lie on those contours which have the correctly quantized 
energy (semiclassical approximation |27j]). In Fig. |^ is 
the situation for potential v; 4> lies on the heavily drawn 
contour shown, which is a bound state oscillating inside 
the well. In Fig. || is the situation for v; the heavy con- 
tour that cj> lies on is of a different nature, not bound in 
any region in space. However its two branches fall away 
from the origin very sharply, so that although they even- 
tually extend out to large x, the projection on the x— 
plane is much denser near the origin than further away. 
This is how such a potential can support a localized den- 
sity near the origin. The phase of the wavefunction in a 
semiclassical view is given by the action integral J p(x)dx 
along the contour and the large step in p that is made 
in a short step in x thus implies that the phase oscillates 
rapidly. The tails of the density, which is the same in 
both cases, arise from fundamentally different processes: 
in the case of the simple sech 2 well (Fig. ^) the tails arise 
from classically forbidden tunneling whereas in the case 
of the divergent potential (Fig. ||) they are classically 
allowed but have exponentially small amplitude. 




described in the text with c = 0.2. The thinner solid line i 
the top figure is the real part of if). 



Finally, we relate the result for v (Eq. (pS})) to that 
obtained from the approach in |2^] (and outlined hi [II Bp . 
Observe that for one electron, the initial conditions on <£> 
required in |^3| are the same as our Eqs. (|4|) and (|^) 
(or, equivalently Eq. (^6|)), each evaluated at t = 0. In 
one-dimension it is then a straightforward exercise to cal- 
culate the terms of Eq. ( p2| ) at t = 0, and, if we disregard 
the boundary condition, we can thus obtain an equation 
for the slope of the potential. This potential gradient 
is consistent with the equation (^8|) evaluated at t = 
that we obtained by the time-dependent Schrodinger's 
equation above: we get here 

Av' = (cV + 2cnh + 2cj + n 2 f) /n 3 (29) 

where n = n(xO), j — j(xO), v — v(xO), c — c(0), and / is 
a constant to be determined by the boundary condition, 
Av — > at oo. However in the present one-electron case 
such a boundary condition cannot be satisfied. 

B. Two non-interacting electrons in one-dimension 

Let 4>i and <p2 represent the initial orbitals for an initial 
state $ and <p\ and 4>2 represent those for another initial 
state $. We choose 

4>i(x) = (j>i{x) expiOi(x) (30) 

where the 9i(x) are real functions. This form guarantees 
the densities are the same initially. The difference of 
initial current between $ and $ is 

Aj =8[(x) \<t>i(x) | 2 + B' 2 {x)\<h[x) | 2 , (31) 
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where the prime indicates differentiation, and so the con- 
dition of equal initial n becomes 

^(e[(x)\M^)\ 2 + 0' 2 (x)\h^)\ 2 ) =0. (32) 
The choice 

9[(x) = c\(f> 2 (x)\ 2 , and ff 2 {x) = -c|&(a)| 2 . (33) 

where c is some constant, ensures eq. |3l] is satisfied. 

To simplify the calculation of the potential gradient in 
Eq. ( |2^ ) further, we take the orbitals <p\ and <p2 to be 
real and take the density to be time-independent. After 
straightforward calculations we arrive at 

Av' = -(?{n'fyl/n + 2&£i(&# + (34) 

This gives the initial gradient of the potential in which $ 
will evolve with the same density as that of $ at t = 0. 

These equation was used to make Fig. |[ for which 
c = 4. The orbitals <fii(x) are just the lowest and first ex- 
cited states of the harmonic oscillator of force constant 
k = 1. Although this kind of potential is not strictly 
allowed because it does not remain finite at oo, we ex- 
pect that the results still hold for well-behaved initial 
states: the difference between the potentials for $ and $ 
vanishes at oo. Moreover, it, or rather, the interacting 
3-D version (Hooke's atom) is instructive for studying 
properties of density functionals (see e.g. [p8| ) because 
an exact solution is known. We can make the dips on the 
right of the figure arbitrarily big, simply by increasing c. 
Note that the alternative orbitals probably do not yield 
an eigenstate of this potential. In the next instant, this 
alternative potential will change, in order to keep the 
density constant. This change can be calculated using 
Van Leeuwen's prescription. 
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